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Section  0.  Introduction 


One  of  the  major  themes  of  research  in  the  area  of  stable  processes 
has  been  the  analysis  of  their  similarities  and  differences  with  Gaussian 
processes,  both  of  which  there  are  many.  A  characteristic  which  stable 
distributions  and  processes  may  possess,  yet  which  their  Gaussian  counter¬ 
parts  may  never  possess,  is  skewness  (i.e.  asymmetry  in  the  distribution). 
Even  though  this  is  the  case,  most  research  in  stable  processes  has  dealt 
only  with  the  symmetric  case.  While  there  is  much  left  to  be  done  in  the 
symmetric  case,  it  behooves  us  to  investigate  the  asymmetric  case  as  well, 
since  it  adds  a  new  dimension  to  the  theory  and  since  the  freedom  to  allow 
skewness  enhances  the  potential  for  application  involving  these  processes. 
This  work  is  a  modest  step  in  that  direction. 

The  outline  of  the  paper  is  as  follows.  In  section  1,  we  define  the 
relevant  terms  and  give  some  necessary  (and  some  unnecessary)  preliminary 
facts.  In  section  2,  we  characterize  all  stable  independent  increment  pro¬ 
cesses,  and  develop  a  Wiener-type  stochastic  integral  with  respect  to  those 
with  location  parameter  zero.  We  also  introduce  a  canonical  independent 
increments  process  with  "maximum  skewness."  We  then,  in  section  3,  represent 
"most"  strictly  stable  processes  as  stochastic  integrals  with  respect  to 
this  canonical  process.  This  theorem  in  some  sense  includes  and  elucidates 
the  known  spectral  representation  theorem  for  symmetric  stable  processes 


Section  1.  Definitions  and  Preliminary  Notions 


-  1 

The  stable  distributions  arise  as  the  only  possible  limiting  distribu¬ 
tions  of  normalized  sums  of  i.i.d.  random  variables.  P.  Levy  discovered 
these  distributions  and  also  computed  their  characteristic  functions.  A 
random  variable  X  is  stable  if  whenever  and  X2  are  independent  copies 
of  X,  and  c^  and  c2  are  any  positive  constants,  there  exist  real  constants 
c>0  and  d  so  that  cX  +  d  has  the  same  distribution  as  c-jX^  +  c2X2.  Equivalent¬ 
ly,  there  should  exist  real  constants  a,  8,  y,  and  v  with  0  <  o  <_  2, 

-1  <_  8  <_  1,  and  v  >  0  so  that  X  has  characteristic  function  $^(t)  =  Eexp(itX) 
of  the  form 

rv|tja(l-ie  tan  sgn(t))  -  iyt  if  a^l 

(1.1)  -log<f>v(t)  =\  7  2 

[v|t|(l+is  ^  sgn(t)  logjt|)  -  iyt  if  a=l 

Here  a  is  called  the  index  of  stability  (in  which  case  we  call  the  distribu¬ 
tion  a-stable),  a  is  a  skewness  parameter,  y  is  a  location  parameter,  and 
v  is  a  scale  parameter.  The  skewness  parameter  6  gives  a  measure  of  how 
much  of  the  Levy-Khintchi ne  jump  measure  is  placed  on  the  positive  and 
negative  half  lines.  For  example,  8  =  +1,  -1,  or  0  according  as  the  measure 
is  concentrated  on  the  right  half-line,  the  left  half-line,  or  is  symmetric. 

For  a  historically  instructive  and  amusing  account  of  these  distributions 
and  their  characteristic  functions,  see  P.  Hall  [7]. 

X  is  called  strictly  stable  if  d  above  is  identically  zero;  or  equivalently, 
if  y  =  0  in  the  case  a  t  1 ,  or  if  8  =  0  in  the  case  a  =  1 .  If  a  >  1 ,  strictly 
a-stable  distributions  have  mean  zero.  Note  that  in  the  case  a  t  1,  all 
stable  distributions  are  just  translations  of  strictly  stable  distributions, 
but  in  the  case  a  =  l,  the  strictly  stable  distributions  can  have  no  skewness  - 
they  are  just  translations  of  symmetric  Cauchy  distributions.  Because  we 
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lose  almost  no  generality  if  we  consider  only  strictly  stable  distributions 
and  processes  in  the  case  a  t  1,  we  will  often  do  this.  However  for  a  =  1, 
we  must  consider  non-strictly  stable  distributions  and  processes  to  maintain 
generality. 

If  X  is  stable  and  satisfies  (1.1)  we  will  write  X  ~  Sa(v, B,u).  Using 
(1.1),  the  reader  may  easily  verify  the  following  two  facts. 

Lemma  1 . 1  If  X  ~  S  (v,6,y)  then  for  real  c,  cX  ~  S  (|c[av,  sgn(c)3,m)  where 

2 

m  *  cu  if  a  f  1,  and  m  =  cy  -  clogjc|  —  ve  if  a  =  1. 


Lemma  1 . 

■1  If 

X-|  and  X2  are 

independent  with  X. 

J 

,  then 

for 

real  c^, 

O 

« 

CVJ 

O 

1 X1  +  C2X2  ' 

,g,p)  where 

(1.2) 

v  = 

|c1!0lv1  +  |c2 

A* 

9 

(1.3) 

3  - 

(c1<a>v1e1  + 

C2<a>v282^v’  and 

(1.4) 

1 

U  =  < 

clyl  +  C2U2 

2 

i  f  a 

t  1 

1 

.clwl  +  C2U2 

L. 

IT 

(v161c1log|c1  j+v 

2e2c2log|c2|) 

if  a 

=  1 

*»  3 

Here,  and  in  the  sequel  we  use  the  convention  that  for  complex  z  =  re 
and  real  p,  z<p>  denotes  rpe19;  so  that  for  x  real,  x<p>  denotes  [x|psgn(x). 

We  shall  study  stable  processes  in  the  remaining  sections.  But  before 
embarking  on  that  course,  we  give  two  propositions  and  an  example,  irrelevant 
to  the  rest  of  the  paper,  included  only  for  whatever  intrinsic  interest 
they  may  have. 

It  is  well-known  that  for  0  <  p  <  1  and  0  <  a  <_  2  that  if  A  is  a  posi¬ 
tive  p-stable  random  variable,  (i.e.  A~S  (v,l,0))  and  X  is  symmetric  a- 
stable  and  independent  of  A,  that  A^aX  is  symmetric  pa-stable  (see  Feller 
[6],  p.  596,  ex.  9).  Here  is  a  slight  generalization. 


Proposition  1.3  Let  A  and  X  be  independent  with  A~S  (v.,1,0)  and 

P  A 

X'Sa(vx,8,0)  where  0  <  p  <  1 ,  0  <  a  <_  2,  and  a  f  1  t  pa.  Then  A1/2X~S  (7,8,0). 

pa  ' 

In  the  case  a  <  1,  if  g  *  -1,  0,  or  1,  then  J  =  8. 

Remark.  The  values  for  the  parameters  v  and  6  are  implicit  in  the  proof,  their 
explicit  mention  being  avoided  for  aesthetic  reasons. 

Proof.  WLOG  assume  that  vA  =  vx  =  1 . 

With  some  manipulation  of  (1.1)  and  application  of  some  standard 

complex  variable  arguments,  one  can  deduce  that  for  complex  w  with  Re  w  >  0, 

(1.5)  E  exp(-wA)  =  exp(-kwp) 

where  k  =  1/cos  and  we  interpret  wp  =  epl°9w  =  |„|Peip*arg(w)  w1th 
|argw|  _<  ir/2.  Letting  w  >  0  in  (1.5)  we  can  deduce  that  A  >  0  a.s.  We 
compute  the  characteristic  function  of  A1/aX.  For  a  t  1 , 

Eexp(itA1/aX)  =  EE(exp(itA1/aX)|A) 

=  E  exp{- | tA1/a|a(l-iB  tan  sgn(t))} 

=  E  exp{.[jt|a(l-i8  tan  2*  sgn(t))]A} 

=  exp{-k[| t|a(l-ie  tan  ^  sgn(t)]p} 

=  exp{-kjt|pa(l+B2tan2  ^)[cos  Pe  -  -,sinPe  sgn(t)]} 

=  exp{-kcosP0(l+82tan2  ^)|t|p0l[l-i  J  tan  sgn(t)]} 

where  6  =  Arctan  (b  tan  -^-)  and  8  =  tanP0/tan(PaTr/2) .  It  is  not  difficult 
to  see  that  [ 6 i  <_  1  in  all  cases  under  consideration.  Note  also  that  if 
a  <  1  and  8  =  ±1 ,  that  0  =  ±  ar,d  8  =  ±1  as  well.  However,  J  t  s  in 
general .  _ 


exists  and  is  given  by 


(1.6)  E|X|p  =  vp/a  c(p)f(p,a,B) 


where 


c(p)  =  2p_1(pf  u‘p'1sin2u  du)-1 

l  n 


,  0  p/2a 

,2i.._2  ait» K  . 


f(p,a,e)  =  r(l-p/a)(l+0^tan^  ^)  cos(^  Arctan(Btan^)) 


Further,  for  fixed  v,  p,  and  a,  E|X|P  is  even  in  0  and  increasing  in  |b|. 

Proof.  First  note  that  X  is  distributed  as  v1/aX0,  where  XQ~Sa(l ,0,0),  and 
hence  that  E[Xjp  =  vp^a  E|XQ|P.  An  application  of  [16,  Thm.  2]  and  a  few 
changes  of  variable  show  that 

00  r 

E | XQ [ p  =  (4 J\fp-1sin2u  du)_12p+1a_1  J0t"p/a‘1  (l-e^cosCst) )dt 
where  we  have  let  s  =  6  tan  Me  now  compute  the  second  integral. 


O  _  00  .eg 

t"p  a"‘(T-e’tcos(st))dt  =  (l-e"tcos(st))(-2.t_p/c‘)  I  -t~  t‘p/ Vt(cos(st)+ssin(st) ) 

)  p  '0  p  *0 

=  0  +  Re  £  f  t‘p/a(l+is)e't^1+1s^dt 
P  Jn 


•  *•  f  l0'TTfj)‘P/V2^ 

=  Re  ^  (l+is)p/a  r(l-p/a) 


r  r(l-p/a)(l+s2)P//2atcos(^Arctan  s) 


This  shows  (1.6).  That  E|X|P  is  even  in  3  is  clear.  To  see  that  it  is 
increasing  in  |3|.  note  first  that  cos  Arc  tan  (3  tan  4p))  = 

cos(^  Arctan(  |e|*  [tan  3p|}),  and  making  the  substitution  e  =  Arctan(  [ e  j  |tar>^() 
we  see  that  e  increases  with  |3|  and  0^6  <  .  Now,  letting 

g ( 0 )  =  f(p,a,3)/r(l-p/a)  =  (l+tan2e)p^2cicos(£  6) 

=  cos(£  0 ) / (cos  e)P//a  =  cos(r6)/cosr0, 

a 

we  can  compute  that  g'(e)  =  r(cos  e)  sin[6(l-r)]  _>  0.  Hence  f(p,  a,  3)  is 
increasing  in  js | ,  as  is  E [X [ p.  ■ 

Example  1 .5  (Chen  &  Shepp  [3])  Let  Z.,  j  *  1,  2,  3,  be  independent,  each 

J 

with  a  S-|  (1,1,0)  distribution.  Set  X  =  Z-j  -  and  Y  =  Z-|  -  Z^.  Clearly 
X  and  Y  are  symmetric  about  0.  Using  Lemma  1.2,  the  reader  may  easily 
verify  that  X  +  Y  is  also  symmetric,  but  not  about  01 
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Section  2.  Independent  Increments  Processes  and  Stochastic  Integrals 

Let  I  be  an  interval  (finite  or  infinite)  contained  in  FR.  A  process 
{ Z( s ) :  sel }  is  said  to  have  independent  increments  if  the  random  variables 
(increments)  Z(t^)-Z(s1),  ....  Z(tn)-Z(sn)  are  independent  whenever 
S1  1  t]  £  s2  <  t^  <  ...  <_  sn  <_  tn  and  the  s^,  t.  are  in  I.  We  characterize 
those  processes  with  stable  increments  as  follows. 


Theorem  2.1  Fix  0  <  a  £  2.  Let  { Z( s ) :  sd}  be  an  independent  increments 

process  with  each  increment  stable  of  index  a.  Then  there  exist  real -valued 

functions,  p,  v,  and  3  on  I  with  v  increasing,  3  measurable,  and  |g|  <_  1 , 

sued  that  each  increment  Z(t)  -  Z(s)  has  scale  parameter  |v(t)-v(s)|,  skew- 
et 


ness  parameter 


6(x)dv(s)/ |v(t)-v(s) i ,  and  location  parameter  p(t)-p(s). 


Conversely,  any  such  choice  of  functions  p,  v,  and  3  defines  such  an  independent 
increments  process  in  the  obvious  way. 


Remarks  (i)  If  v(t)  -  v(s)  =  0,  the  skewness  parameter  above  is  not  defined, 
but  is  immaterial  -  we  may  take  it  to  be  zero. 

(ii)  This  theorem  has  an  obvious  analog  in  the  setting  of  independently 
scattered  stable  measures. 


Proof.  For  any  stable  variable  X,  let  us  define  V ( X )  =  v,  B( X)  =  3  and 
M(X)  =  p,  where  v,  3,  p  are  as  in  (1.1).  Now  pick  tQcI  and  set  Z'(t)  = 

Z(t)  -  Z(tg)  for  all  tel.  Since  each  Z'(t)  and  each  increment  Z'(t)  -  Z'(s) 
is  stable,  we  may  define  the  functions  v(t)  =  V ( Z 1 (t) )sgn(t-tc),  b(t)  = 

B ( Z ' (t) ) ,  w(t)  =  b(t)v(t) ,  and  p(t)  =  M(Z'(t)). 

We  claim  that  for  s  <_  t,  v(t)  -  v(s)  =  V( Z(t)-Z(s ) ) ,  and  thus  v  is 
increasing  and  Z(t)  -  Z(s)  has  scale  parameter  jv( t)-v(s) j  for  arbitrary 
s,t.  We  verify  the  claim  only  in  the  case  s  £  t  <_  tg,  leaving  the  similarly 
proved  cases  s  _<  tQ  ^  t  and  <  s  <  t  to  the  reader.  In  the  case  at  hand. 


V(t)  -  v(s)  =  V(Z(t)-Z(t0))sgn(t-t0)  -  V( Z(s )-Z( tQ) )sgn(s-tQ) 

=  V(Z(s)-Z(t0))  -  V(Z(t)-Z(t0) ) 

=  V(Z(t0)  -  Z(t)  +  Z(t)  -  Z(s))  -  V(Z(t0)-Z(t)) 

=  [V(Z(t0)-Z(t))+V(Z(t)-Z(s))3  -  V(Z(t0)-Z(t))  by  (1.2) 
=  V(Z(t)-Z(s) ) 

Since  v  is  increasing,  it  is  measurable,  and  induces  a  Lebesgue-Stieltjes 
measure  dv  on  I. 

Now  note  that  by  Lemma  1.2,  if  and  X2  are  independent  and  stable, 

that 

b(x1+x2)v(x1+x2)  =  b(x1 )v(x1 )  +  b(x2)v(x2) 

Hence  for  tQ  _<  s  <_  t,  letting  X1  =  Z(t)  -  Z(s)  and  X2  =  Z(s)  -  Z(tQ), 
w ( t )  -  w(s)  =  b(t)v(t)  -  b(s)v(s) 

=  b(x1+x2)v(x1+x2)  -  B(X2)V(X2) 

*  B(X1)V(X1) 

-  B(Z(t)-Z(s))[v(t)-v(s)] 

It  is  easy  to  verify  in  like  fashion  that  the  equation 
w(t)  -  w(s)  =  B(Z(t)-Z(s))[v(t)-v(s)] 
holds  true  in  the  cases  s  _<  tQ  _<  t  and  s  £  t  _<  tQ.  Hence 

(2.1)  w(t)  -  w(s)  =  B(Z(t)-Z(s))!v(t)-v(s)| 


for  arbitrary  s,t. 

In  any  case,  |w(t)-w(s)  j  _<  v(t)-v(s)  for  s  _<  t.  Hence  for 

n  n 

a  1  1  ti  **•  i  tn  i  b»  }  iw(tj)-w(tj_1)|  i  z  v(tj)  -  v( t^^ )  <_  v(b)  -  v(a). 

J  1  0  1 

This  shows  that  w  is  of  bounded  variation  on  any  finite  interval  contained 
in  I  and  induces  a  Lebesgue  -  Stieltjes  measure  dw  on  I  which  is  absolutely 
continuous  with  respect  to  dv.  We  may  set  s(x)  =  ^  (x).  Clearly  |g(x)|  <_  1 


a.e.  [dv]  and  by  (2.1),  j  B(x)dv(x)/ jv(t)-v(s) j  =  B(Z(t)-Z(s) ) ,  which  is 

Js 

the  skewness  parameter  of  Z(t)  -  Z(s). 

Finally,  it  is  easy  to  check  (using  Lemma  1.2)  that  M(Z(t)-Z(s))  = 
u( t)  -  y(s)  for  all  s,t.  This  holds  even  in  the  case  a  =  1,  since  log | ±1 j  =  0. 

We  now  introduce  a  particular  independent  increments  process  of  some 
importance  to  us  in  Section  3.  Fix  a  and  let  Va(t)  be  the  right-hand  side 
of  (1.1)  with  v  =  3  =  1  and  p  =  0.  Since  exp(-i^(t))  is  the  characteristic 
function  of  an  infinitely  divisible  distribution,  there  is  an  independent 
increments  process  {R(s ) :  sefRlwith  characteristic  function 

E  exp(itR(s) )  =  exp(-|s|i|>o(t)) 

(see,  e.g.[l]).  This  process  R  will  be  called  the  canonical  totally  right- 
skewed  a-stable  Levy  process.  The  corresponding  function  parameters  of 
Theorem  2.1  for  R  are  v(t)  =  t,  6( t )  a  0.  For  0  <  a  <  1 ,  R  is  an  increasing 
process. 

We  now  wish  to  produce  a  Wiener-type  stochastic  integral  with  respect 

a  process  with  independent  stable  increments.  The  development  here  is  standard, 

but  some  interesting  things  happen  in  the  case  a  =  1. 

To  lay  the  groundwork  for  the  integral  in  this  case  we  introduce  the 

Orlicz  space  L  log+L,  defined  as  follows.  Let  (M,dp)  be  a  finite  measure 

space  and  let  $(x)  =  x  log+x  for  x  _>  0.  Note  that  $  is  positive,  convex, 

and  increasing  on  [0,®),  and  that  $  satisfies  the  so-called  condition  at 

°°,  i.e.  lim  <t>(2x)/<t(x)  <  «.  This  is  enough  to  guarantee  that  the  Orlicz 
x-*» 

space  based  on  4>,  which  we  call  L  log+L(M,dp),  defined  by 

L  log+L  =  {measurable  f:M-*fR:  i  $(  jf j )dp<~}, 

JM 


is  a  Banach  space  under  the  norm 


li^llLlog+L  =  inf(c>0:  j  <f(  |f  | /c )d y<1 } 

M 

and  that,  in  this  space,  simple  functions  form  a  dense  subset. 

What  is  important  for  us  is  that  convergence  of,  say,  f  to  f  in  Llog+L 


implies  the  same  convergence  in  L  ,  and  that  of  the  integrals  j ( fp-f )log j f n-f I 
to  zero,  as  is  shown  in  the  following 


Lemma  2.2  For  Ilf  II, ,  +,  =  c  <  e' 

-  Llog  L  — 


(2.2)  II f II 1  <_  K,c 


(2.3)  J f  1  og  |f  j  jd m  <_  K7Jc  |logc  J 


the  constants  K-j  and  l<2  depending  only  on  u(M), 


Proof  Note  that  by  the  monotone  convergence  theorem,  we  have  that 


ill  l0g+(M)  <1.  Hence 


1  >  f  ifl  !0,  ill  >  f  iii 

—  c  3  c  —  c 


so  that 


{ |  f  |  >c } 

(  ,  , 

I f I  ic. 

{ |  f 1 >ce} 


{ If  i >ce) 


Since  J  | f j  _<  cey(M),  (2.2)  is  proved  with  =  1  +  eH(M). 

{|f |<ce) 

For  (2.3),  examination  of  |xlogxj  for  xc[0,l]  shows  that  for  6<  e'i 


(2.4) 


flog|f||  = 


^f  ll 


{  |f  |  >1  }  {frjfljil}  {  ]f  i<6} 


|  f |1 og  j  f  j  +  y{ j  f | >6 }e” '  +  6|log6ju(M) 


which  shows  that  j  |fj1og|f|  ±c.  Noting  that  y{|f|>6}  <  llfll^-1  by  the 

t|f|>i> 

Markov  inequality,  and  choosing  6  =  ^c  gives  us  through  (2.2)  and  (2.4)  that 

|  |f  log|f  1 1  <_  c  +  ilfj  U6”^e~*  +  ^(Mj^c  1 1  og  c  j 

1  c  +  Sc  (e_1+u(M))  +  j  u(M)v^c|log  c|, 

establishing  (2.3).  ■ 

Now,  let  (Z(s):  sel)  be  a  process  with  a-stable  independent  increments 
as  in  Theorem  2.1.  For  ease  and  simplicity  we  will  let  Z  be  such  that 
u(t)  s  0,  although  processes  with  sufficiently  well-behaved  y(t)  can  be 
handled  in  a  straightforward  way.  Also,  we  assume  in  the  case  a  =  1  that 
v (I )  <  °°. 

n 

For  a  step  function  f  of  the  form  f  =  z  c.l/  (where  c.dS, 

j=l  J  lSj-l)sjJ  J 

sjd,  sj_i<Sj  for  all  j)  we  define 

(2.5)  S(f)  =  f (s)dZ(s)  =  l  c.[Z(Si)-Z(s.  ,)] 

Jl  j=1  J  J  J-i 

Theorem  2.3  The  stochastic  integral  S(f)  =  f(s)dZ(s)  can  be  defined  for 

-*1 

all  f  in  the  space  La(I,dv)  (resp.  Llog  +L(I,dv))  in  the  case  a^l  (resp. 
a=l)  so  that  it  is  linear,  agrees  with  (2.5)  for  step  functions,  and 
satisfies  S(f)'Sa(v^,6f,yf)  where 


vf = 


h 


|f(x)|adv(x)  =  llfll*  , 


■  1 


[f(x)]<a>S(x)dv(x), 


f(x)logjf(x) |e(x)dv(x) 


and 

tfl 

a=l 


Further  $(f)  and  S(g)  are  independent  if  and  only  if  fg  =  0  a.e.  [dv]. 


Proof  Note  that  by  Lemma  1.2  and  Theorem  2.1,  for  a  step  function  f  as  in 
(2.5)  we  have  S(f )~Sa(vf, 6f,yf)  where 

»f  *  £  lcj l“Cv(tJ)-v(tj_1 )3 

J  * 

t  j 


Sf  =  Vf"1  E  cj<a>[v(tj)-v(tj_i)]  eCxJdvCxJ/tvCt^-vCtj  .,)] 

j*l 

*  i  1 


J 

'j-1 


v;1  f  f<a>e  dv 


and 


a/*l 


uf  - 


U  j=lCjl°9lcj  iCv(ti  )-v(tj_1 )]  3(x)dv(x)/[v(t1  )-v(tj._1 )]  a=l 

Vl 


f  log |f |gdv 


a>0 

a=l 


I 


Thus  the  claim  concerning  the  distribution  of  S(f)  is  true  for  step  functions. 
Also,  for  such  functions,  the  claim  of  linearity  is  easily  seen  to  be  true. 

To  define  S(f)  for  general  f  in  La  (resp  L  log+L),  pick  a  sequence  of 
step  functions  f  converging  in  norm  (quasi-norm,  if  a<l )  to  f.  Note  that 


r  m  » 
r'- 


V*  i’'  /’**'  »*■  v  ••  vV  ' 


■ 

f  dZ  -  f  dZ  =  S(f  -f_)  by  linearity  on  step  functions,  and  so 
n  m  n  ni  '  r 

f  a 

fRdZ  -  fmdZ  has  scale  parameter  v^  =  II f n- fm II a  and  location  parameter 

J  J  n  m 

2 

Mf  s  0  if  a  i1  1  and  “  -  —  ( ^  - ^  )  1  og j -f  - f  1  edv  if  a  =  1.  Since  these 

n  m  n  m  J  "  111  n  m 

parameters  approach  zero  as  n,m  become  large  (cf.  Lemma  2.2  in  the  case 

a=l ) ,  we  conclude  that  {jf  dZ}  is  a  Cauchy  sequence  in  the  metric  of  con- 

n  r  r 

vergence  in  probability.  We  are  then  free  to  define  i  fdZ  =  lim  prob  j  f  dZ. 

j  n-*=°  n 

Since  Ilf  II  -  Ilf  II  ,  ['  f  <a>3dv  -*•  f  f<a>3  dv,  and  f  f  log  If  igdv  -*• 
na  a  J  n  J  j  n  * 1  rv 

f  1 og  j f 1 6dv ,  Levy's  continuity  theorem  guarantees  that  and 

are  as  advertised.  S(*)  is  linear  by  continuity. 

To  verify  the  claim  concerning  independence,  we  observe  that  S(f)  and 
S(g)  are  independent  if  and  only  if 


*S(f),S(g)(s,t)  =  *S(f)(s)'*S(g)(t) 

for  all  s  and  t.  If  fg  =  0  a.e.  [dv],  one  can  verify  easily,  using  the 
form  of  the  characteristic  function  for  these  stochastic  integrals,  that 
this  equation  holds.  Conversely,  for  S(f)  to  be  independent  of  S(g),  we 
must  have  that 


l*S(f),S(g)(s,t)l  =  '9S(f)(s)il*S(g)(t)i 
which  implies  that 


lisf+tglla  =  llsf#“  +  ll  tg  II a 


Section  3.  Stable  Processes  and  their  Spectral  Representations 


A  stochastic  process  iX^itcT;  will  be  callec  a  s^aule  ^rosess  if  each 

random  vector  Y  =  (X.  ),  t.<.T,  is  multivariate  stable,  i.e.  if  for 

_1  _  n  _  J 

independent  copies  X^  and  of  X  and  positive  constants  c-|  and  C2,  there 

exist  constants  c  >  0  and  dt.Wn  so  that  c-jY-,  +  0^2  is  distributed  as  cY+d. 

As  before,  tX  }  is  a  strijilj  stable  process  if  d"  is  zero  for  all  t.  and  c. 

v  J  J 

Clearly  if  Y  is  stable,  then  *-Y  is  stable  for  any  7<.fRn.  The  converse  is 
not  true  in  general  (see  [13]),  but  is  true  in  the  strictly  stable  case 
(see  [4]). 

Levy  [12]  and  Feldheim  [5]  have  shown  that  a  random  vector  Y  in  Rln  is 
multivariate  stable  if  and  only  if  there  exists  cu.(0,2],  peFRn,  and  a 
finite  measure  r  on  the  unit  sphere  Sn  of  IKn  such  that  the  characteristic 


function  ^  of  X  satisfies 


(3.1)  -log  *y(t) 


Is"1™*- 

1  (it -sj  + 


i  tan  -y-  (t*s)<a>dr(s)-it*u  if  a^l 

|-(‘t*?)log jt*T| )dr(T)  -  it*7  if  a=l 


When  n  =  1,  this  agrees  with  (1.1)  with  v  =  '.'(S1  )  and  z  =  [;  (1 )-?(-!  )]/?(S^ ) 
Also  observe  that  this  shows  that  the  random  variable  t-Y  has  a  S  (v, ,>,•„) 
distribution,  where  v  =  I i t-?jadr(?),  ,  =  V_1  ’(t-s) <u>d;(?) ,  and  ,  =  t-~  in 


case  a/1  and  p  =  t*p  -  |  ‘(t-sjlog^s^i'fs)  if 


We  warn  the  reader  that  in  case  a  =  1  this  shows  7  of  (3.1)  is  not  the  vector 

of  location  parameters  of  the  individual  X.'s,  but  rather  X.  has  location 

2  f  J  J 

parameter  p -  r  s  logjs, ,dr(s). 


In  the  representation  (3.1),  strict  stability  is  equivalent  to 
_  ’ _ _  _ 

y  *  0  if  a  /  1,  and  to  (t*s)dr(s)  ■  0  for  all  t  if  a  =  1.  The  latter 

* 

occurs  if  and  only  if  we  may  choose  r  symmetric.  Also,  X  is  symmetric 
about  0  (for  any  a  <  2)  if  any  only  if  IT  =  0  and  we  may  choose  r  symmetric 
Now,  for  an  arbitrary  index  set  T,  pick  arbitrary  functions  {f^iteT} 
in  La([0,l],x)  if  a  f  1  (resp.  Llog+L([0,l],x)  if  a  =  1)  and  set 


Xt  =  f  ft(s)dR(s) 
[0,1] 


where  R  is  the  canonical  totally  right-skewed  a-stable  Levy  process  (see 
section  2).  We  now  argue  that  { >  is  an  a-stable  process  (strict,  if  a^l ) 
Let  R  be  the  canonical  right-skewed  a-stable  Levy  process  on  [0,2].  Let 


ft,i(x)  ■ 


ft(x)  xe[0,l ] 


and 


xe(l,2] 

xe[0,l] 


ft,2(x)  =  \ 


tft(x-l)  xe(l ,2] . 


Cl 


early  the  processes  (Xt  and  {X^  defined  by  Xt  1  =  Jft  and 


t,2 


ft  2^R  are  independent  and  distributed  as  { Xt }  by  Theorem  2.3. 


Straightforward  (and  arduous,  if  a*l )  computations  with  the  joint  char¬ 


acteristic  function  and  formulae  from  Theorem  2.3  show  for  c-j,  >  0 


hat  c  ^  ( X^  i  $  •  •  • »  X^  -j)  +  c2^xt  2  *  *  ■  *  X^  2^  isdi 

c1a+c2a)1^(Xt  .....  Xt  )  +  d ,  where  d”  =  0  when  a  f  1 , 


that 

( 


stributed  as 
and  d^  = 


f  ('  2 

k(ift  dx,  ...  ft  dx)  with  k  *  -[(c1+c2)log(c1+c2)-c1logc1-c2logc2]  when 
*  n  un 

a*l .  Hence,  { Xt }  is  in  fact  a  stable  process.  It  is  always  strictly 
stable  if  a^l  and  is  strictly  stable  in  the  case  a*l  if  and  only  if 


jf^dA  =  0  for  all  t. 
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It  is  the  purpose  of  this  section  to  show  that  practically  all  stable 
processes  can  be  represented  in  this  way,  i.e.  as  a  set  of  stochastic  integrals 
of  appropriate  functions  against  this  canonical  totally  right-skewed  te¬ 
stable  Levy  process.  Our  method  is  in  essence  the  same  as  that  of  Kanter 

in  the  case  of  symmetric  stable  processes  (see  [9]  in  view  of  [14]).  That 

is,  we  shall  use  the  finite-dimensional  representation  (3.1)  in  conjunction 
with  an  imbedding  theorem  (Theorem  3.1)  similar  to  the  main  theorem  of  [9], 

Me  comment  that  another  route  to  this  representation  is  to  use  the  Hilbert 

space  analogue  of  equation  (3.1)  due  to  Keulbs  [10]  and  to  proceed  as  he 

has  in  his  last  section,  carrying  the  skewness  parameters  along.  We  prefer 
to  proceed  in  the  manner  described  because  we  feel  it  is  more  basic  and  direct, 
and  allows  us  to  present  Theorem  3.1. 

Theorem  3,1  (i)  Let  0  <  p  <  °°.  Let  L  be  a  (real  or  complex)  separable 

Frechet  space  such  that  for  each  finite-dimensional  subspace  M  there  is  a 
linear  isometry  JM:  M  -*•  LP[0,1].  Then  all  of  L  is  linearly  isometrical ly 
imbeddable  in  Lp  via  J:  L  ■*  LP[0,1]. 

f 

(ii)  If  in  addition  we  assume  that  for  fixed  f,  !(JMf)<p>dA  is  inde¬ 
pendent  of  M  containing  f,  we  can  then  choose  J  above  so  that 


(3.2)  |( Jf )<p>dA  =  j(JMf)<p>dA 

for  all  f  in  L  and  all  M  containing  f. 


Remarks  Part  (i)  of  the  theorem  in  the  real  case  is  due  to  Bretagnolle 
et  al  [2]  in  the  case  p  _>  1  and  to  Schreiber  [15]  in  the  case  0  <  p  <  1 . 


Both  these  papers  prove  this  without  the  separability  hypothesis  on  L 
and  with  the  range  of  J  being  Lp(E,u)  with  ( E , ^)  a  suitable  measure  space. 
Kanter  [9]  has  given  a  shorter  proof  in  the  real  separable  case  for  0  <  p  £  2 
Part  (ii)  is  our  main  concern  here  -  part  (i)  is  included  because  it  follows 
easily  from  our  proof  (by  ignoring  the  extra  assumption  in  (ii)  and  its 
consequences),  and  because  our  proof  is  a  still  simpler  and  more  general 
version  of  Kanter's. 


Proof  Since  L  is  separable,  we  may  find  f^.fg,  ...  in  L  such  that 
sp{f,,f2,...}  =  L  and  z tl f n»  <_  1  (where  IMI  is  the  quasi-norm  on  L). 


Define 


Mn  =  sp{f-j,...» f n }  and  let  Jn  be  the  hypothesized  JM  . 


Setting  r  =  pvl  we  have  that  for  n  _>  m  and  scalars  A., 

J 


(3.3) 


I  I  A  .  J  f.llP 
3=1  J  n  J  P 


i,W 


0*1 


m 


■  2  A  .  J  f  . 


U  A  .  \J  .  II  ft 

j_l  J  m  J  p 


and  also  (under  the  assumption  of  (ii)) 


(3.4) 


m  -  m 

‘  “oVo1  ‘  j(1LV«fj)  di 

J  '  w 


Call  Fn  =  (  x  |Jnf.jp)1/p  and  define  dyfl  =  Fpd a.  By  (3.3)  and  (3.4), 


j=l 


n  j 


(3.5) 
and 

(3.6) 


m 


r  m  r  - 

'  z  A.J  f./F  jpdy  =  ;  Z  A . J  f ./F  i pdu 

j.j  J  n  j'  n1  Mn  ;  'j=1  j  m  j'  m1  m 


•  m  r  m 

J  I  J  1 


Now  call  B  =  {(z1  ,z2, . .  ,)eK*:i|Zj  |P_<1 ;  where  K  is  the  scalar  field  ® 
or  CC.and  give  B  the  relative  topology  as  a  subspace  of  the  product  space 


D  =  n  {zeK :  jz  |_<1 }.  Since  B  is  closed  in  D,  and  0  is  compact  (by  the  Tychonoff 
j=l 

theorem),  we  conclude  B  is  compact. 


Define  measures  on  B  by 
n  J 


VE>  -  "n{<WF 

for  Borel  sets  E  of  B.  With  this  notation,  equations  (3.5)  and  (3.6)  now 
read 


(3.7) 


m 


I  E  A.z.|pdv „(Z)  =  i  z  X.z.|pdv  (z) 
j=l  J  J  n  JB  j=l  J  J  m 


m 


and 


(3.8) 


m  r  nn 

(  z  X.z.)<p>dvn(z)  =  (  Z  x.z  )<p>dvm(z) 

B  j-1  J  J  n  JB  j-1  J  J  m 


Note  that  vn(B)  =  un([0,l]) 


f1  d  n  D  n  r 

FpdX  *  Z  IIJJ\.IIP  =  Z  llfjlr  <  1.  Since 


0  "  j-1  n  J  p  j=l  J  ~ 

B  is  compact,  the  Banach-Alaoglu  theorem  guarantees  the  existence  of  a 


subsequence  {v  }  and  a  finite  measure  v  on  B  with  {v  }  converging  weakly 
"k  m  mnk 


to  v.  Since  the  functions  z  •»  [  z  X.z.|p  and  z  -+  (  z  X.z.)<P>  are  continuous 

j-1  J  J  j-1  J  J 

on  B  for  each  m  and  X,,...,  x  cK,  it  follows  from  (3.7)  and  (3.8)  that 

I  III  \ 


(3.9) 


J 


m  m  r 

z  x.z.  |pdv(z)  =  II  z  x.f.li 


i-1  J  J 


J  J 


and 


(3.10) 


m 


m 


(  z  x.z.)<p>dv(z)  *  f(j  (  z  x.f.))<p>dx 

j-1  J  J  J  M  j=l  J  J 


for  any  M  containing  M  . 


Now  define  U:  sp{f1 ,f2, 


LK(B,dv)  to  be  the  linear  extension  of 


the  map  fn  pn  where  pn(z)  =  z^.  By  (3.9),  U  is  an  isometry  and  hence  has 
a  unique  extension  (also  called  U)  by  continuity  to  all  of  L.  Also,  under 
the  assumption  of  (ii),  this  extension  must  have  the  property  that 


(3.11) 


|(Uf)<p>dv  =  j(JMf)<P>dx 


for  all  f  and  finite-dimensional  M  containing  f,  by  (3.10).  By  using  a  measure 
algebra  isomorphism  of  B  into  [0,1],  we  may  find  an  isometry  V  of  Lp(B,dv) 

f 

into  Lp( [0,1  ] ,  A)  which  satisfies  j ( Vf ) <p>dA  =  f<p>dv.  Setting  J  =  VU 
completes  the  proof.  ■ 

With  this  theorem  in  hand  we  can  prove  the  representation  theorem 
previously  advertised. 

Let  {X^:teT}  be  an  arbitrary  stochastic  process.  We  say  that  it  satisfies 
condition  S  if  there  exists  a  countable  set  T'  £  T  such  that  every  Xt  is 
a  limit  in  probability  of  a  sequence  from  the  set  of  all  finite  linear  com¬ 
binations  EA  .X.  ;  A.eIR  ,  t.eT'  (see  [8]  and  [10]).  Practically  all  processes 
J  j  J  J 

of  interest  will  satisfy  this  condition,  including  those  continuous  in 
probability  and  those  defined  on  a  countable  index  set  T. 

Let  L^  be  the  linear  span  of  the  X^'s,  closed  under  the  metric  of  con¬ 
vergence  in  probability.  It  is  not  a  difficult  exercise  to  show  that  condition 
S  is  equivalent  to  the  requirement  that  L ^  be  separable. 


Theorem  3.2  Let  a  t  1.  {Xt:  teT}  is  a  strictly  stable  process  of  index  a 
satisfying  condition  S  if  and  only  if  there  exists  a  collection  of  functions 


{ft:  teT}  £  Lot[0, 1  ]  such  that  the  process  {Y^}  defined  by  Y^  =  f^dR  is 


distributionally  equivalent  to  {X*}. 


Proof  For  the  "if"  part  of  the  theorem,  we  note  that  {Yt>  above  is  a  strictly 
oi-stable  process  by  arguments  given  at  the  beginning  of  this  section.  To 
see  that  {Y^}  (and  hence  {X^ >)  satisfies  condition  S,  note  that  La[0,l] 

(and  hence  any  subspace  thereof)  is  separable,  and  the  closed  linear  extension 
of  the  map  ft  -*•  Y^  is  an  isomorphism  of  spif^a  onto  Ly. 

We  now  prove  the  "only  if"  part.  Clearly,  every  element  of  I_x  is 
strictly  stable  of  index  a,  and  moreover,  {ZiZeL^}  is  a  strictly  a-stable 
process.  Now  let  M  be  a  finite-dimensional  subspace  of  L x,  and  let  Z^,...  ,  Z 
be  such  that  M  =  sp{Zj,...,  Z^}.  Since  I  =  (Z, ,  ....  Z  )  is  jointly  strictly 
a-stable,  we  may  write  by  (3.1). 


(3.12) 


-log^t)  =  ( i t -s  |a-i  tan  ^  (t*s)<ci>)dr(s) . 


for  some  finite  measure  r  on  S 


n 


Define  g.:  Sn  -*•  JR  by  g.(s)  »  s  .(r(Sn))'1^“.  Now  pick  a  measure-pre- 

J  J  J 

serving  set  transformation  T  of  (F,(r(Sn))’1dr)  into  (8, d A ) .  Here, 

F  =  :  j  -  1.  ....  n},  8  =  {Borel  sets  of  [0,1]}, 

and  dx  is  Lebesgue  measure.  Set  h.  =  Tg..  It  is  evident  from  (3.12)  that 

J  J 


(3.13)  -log^ft)  =  j  ( iztj9j(s)  |a-i  tan  f  (z  tj9j (?) )<a>) (r(Sn) )*1dr(?) 

Sn 


A 


r« 


j  jzt.h .(x) jadx(x)  -  i  tan  —  |  (it .h .(x))va>dx(x) 
JfjJJ  ^  Jq  J  J 


Now  we  metrize  Lx  as  follows.  For  any  YeL^  we  set  IIYil  =  (-log |EeiY | J1 
By  (3.13),  IMI  a  quasi-norm  (a  true  norm,  if  a>l )  on  L x  which  metrizes  con¬ 
vergence  in  probability.  Also,  for  so  normed,  the  linear  extension  of 
the  map  Z.  -*•  h.,  call  it  JM,  is  an  isometry  of  M  into  La[0,l].  We  also  see 

j  j  i  «  < 

that  for  each  Z  in  M,  j  (0M(Z))<a  dx  is  just  the  skewness  parameter  of  Z, 

J  o  n 


dv1 

multiplied  by  liZIl  ,  and  hence  does  not  depend  on  M  containing  f. 

We  can  now  apply  Theorem  3.1  (ii)  to  conclude  that  there  is  an  isometry 
J:  U-H-a[0,1]  with  (J(Z))<a>dX  being  the  skewness  parameter  of  Z,  multiplied 

*  i 

by  HZIiav  .  Setting  =  J(Xt)  and  recalling  Theorem  2.3  completes  the  proof.  ■ 

The  author  has  not  as  of  yet  been  able  to  obtain  an  analogous  version 

of  this  theorem  in  the  case  a  =  1.  One  difficulty  stems  in  part  from  the 

curious  behavior  of  the  location  parameter  (cf.  Lemma  1.1,  (1.4),  example  1.5, 

Theorem  2.3,  and  the  comments  following  (3.1)).  It  is  not  a  linear  function 

of  the  variables  in  the  process  as  is  the  case  if  a  M.  Another  difficulty 

comes  from  the  fact  that  with  respect  to  scale  and  skewness,  the  appropriate 

"imbedding  space"  is  L^;  yet  for  the  location  parameter  to  exist,  integrals 

of  the  form  f  log | f  jd A.  must  be  defined  in  some  appropriate  sense.  Hence 

. 

it  is  not  clear  what  the  appropriate  space  of  integrable  functions  should 

be  -  L1  is  too  large  (but  has  the  "right"  norm  in  one  sense),  and  Llog+L 

may  be  too  small  or  may  have  an  altogether  inappropriate  norm. 

We  close  by  making  a  few  comments  regarding  Theorem  3.2.  First,  we 

observe  we  have  lost  little  generality  in  assuming  strict  stability,  since 

subtracting  location  parameters  (which  are  means  if  a>l )  from  a  general 

stable  process  yields  a  strictly  stable  process  in  the  case  a^l .  Secondly, 

we  mention  that  in  some  sense  the  skewness  of  elements  of  the  process  are 

mirrored  in  the  skewness  of  the  corresponding  representing  functions.  More 

specifically,  observe  that  a  process  Xt  =  j  ftdR  has  all  variables  totally 

right-skewed  (resp.  left-skewed),  i.e.,  all  variables  have  skewness  parameter 

+1  (resp.  -1),  if  and  only  if  f^  _>  0  (resp.  f^O)  a.e.  [dx]  for  each  t.  Also  if 

{X. }  is  a  symmetric  stable  process,  setting  r  ■> 

z  „  \  ft(2x)  for  0  <_x  <  j 
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and  observing  that  { j f t d R >  is  distributed  as  {X^}  gives  us  that  iX^.}  is  a 
symmetric  stable  process  if  and  only  if  it  has  a  representation  { j f ^dR } 
with  each  ft  (and  each  linear  combination  ZA.f.  )  symmetrically  distributed. 
Also,  our  theorem  relates  to  the  spectral  representation  theorem  for  symmetric 
stable  processes  (see  [8])  as  follows:  if  the  symmetric  a-stable  process 


iX^}  is  represented  by  { j f ^dZ > ,  where  Z  is  a  symmetric  a-stable  Levy  process, 
then  defining  f^  from  f^  as  above,  the  process  {jftdR>  also  represents  L Xt > . 
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